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Finite-temperature spin dynamics in planar t — J model is studied using the 
method based on the Lanczos diagonalization of small systems. Dynamical spin 
structure factor at moderate dopings shows the coexistence of free-fermion-like and 
spin-fluctuation timescales. At T < /, the low-frequency and static susceptibil- 
ity show pronounced T dependence, supporting a scenario, related to the marginal 
Fermi-liquid one, for the explanation of neutron-scattering and NMR-relaxation ex- 
periments in cuprates. Calculated NMR relaxation rates reasonably reproduce ex- 
perimental ones. 
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The understanding of spin and charge dynamics in strongly correlated systems, as realized 
in cuprates on doping the reference antiferromagnetic (AFM) insulator, still represents a 
challenge for theoreticians. The low-frequency spin dynamics and static spin response in the 
undoped and doped AFM state in cuprates have been recently extensively studied by neutron 
scattering [1] and in a number of NMR and NQR experiments [2,3]. They established that 
normal-state spin dynamics differs qualitatively from the one expected for Landau Fermi 
liquids. Generically, NMR and NQR spin-lattice relaxation time Ti is nearly T and doping 
independent in the normal state T > [3], in contrast to the Korringa law T^^ oc T 
in normal metals. Also, in the same regime the low-cu dynamical susceptibility in doped 
systems appears to be consistent with x"(^) ^ io/T [1]. 

Well understood so far are only undoped cuprates, which behave in all respects as 
isotropic quantum AFM, with the long range order at T = 0. For doped systems NMR 
and NQR data on spin dynamics have been interpreted within the phenomenological model 
of AFM correlated spins, where the T dependence is attributed to the variation of AFM 
correlation length ^{T) [4]. At low hole doping dynamics has been mapped to related 
quantum-critical (QC) scaling regime of nonlinear sigma model where ^ oc 1/T [5] . An 
alternative scenario for the low-cu, low-T behavior has been given in terms of an anomalous 
T-dependence (not directly related to (^(T)), introduced within the marginal Fermi-liquid 
hypothesis [6]. 

Among the microscopical models the most helpful for the discussion of spin dynamics in 
strongly correlated metals is the t — J model [7] 



where c-,.{cis) are projected fermionic operators, taking into account that the double oc- 
cupancy of sites is not allowed. In spite of its simple form the model proved to be very 
difficult to analyse, both analytically [7] as well as numerically [8]. Most reliable results 
related to spin dynamics have been obtained by exact diagonalization studies [8-10] of small 
systems, and by means of the high-temperature series expansion [11,12]. These methods. 
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however, have not been able so far to get results for low-cu and d.c. spin response in the 
most challenging parameter regime J < t and at intermediate < T < J . 

Recently present authors introduced a new numerical method, based on the Lanczos 
diagonalization of small systems combined with the random sampling over basis states [13], 
which allows the study of dynamical and d.c. response functions at T > 0. The method has 
been already applied to the evaluation of optical and d.c. conductivity in the t — J model 
on a square lattice [14]. 

In this paper we present results, obtained by the same method, for the dynamical spin- 
structure factor (we furtheron choose units with h = ks = 1) 



Since numerical requirements are the same as for the conductivity problem, we refer to the 
description of technical details to Ref. [13,14]. We study the planar system of = 4 X 4 = 16 
sites with J/t = 0.3 and variable doping, i.e. with = — 10 holes. Typically we use 
up to Af = 120 Lanczos steps and the sampling over A^o ~ 300 states. The method has 
been tested with the fuU-diagonalization results for = 10 [9]. In comparison with latter 
results our discrete spectra obtained for = 16 are much more dense, so that even minimal 
additional broadening (with a characteristic width rj <^ J) yields smooth macroscopic-like 
spectra. Generally, results become unrealistic at low temperatures T < T* , where finite-size 
effects start to introduce size-dependent features. T* is related to the average level distance 
in the low-energy sector. The latter is smallest (for fixed A^) in the intermediate-doping 
regime 0.12 < < 0.5, where we reach T* ~ O.lt. The low-energy sector becomes more 
sparse both in the low-electron > 0.5 as well as in the pure AFM regime < 0.12, 
leading to the increase of T*. The onset of finite-size effects is monitored by the appearance 
of unphysical structures in the high-frequency spectra, by the dependence on smoothing rj 
etc. We choose furtheron predominantly rj = 0.07 t. 





and related dynamical susceptibility x((f, cu) 




(3) 
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We present in Fig. 1 S{Q^iLi) spectra at the AFM wavector Q = (7r,7r) and fixed T = 
0.2 t < J , but varying hole concentration c/^. As noted above, moderate doping results are 
reliable in this T < J regime, while low-doping spectra already exhibit finite-size effects, 
e.g. gaps appearing in spectra etc. The most interesting feature in Fig. f is the qualitative 
change of spectra on doping. While at < 0.f2 spectra are dominated by a single central 
peak with the width ~ 2J due to AFM fiuctuations, in the intermediate regime 0.f2 < 
c/i < 0.5 a high-frequency component with ~ t emerges, coexisting with the remaining 
low-frequency fluctuations. It is plausible to attribute the high-cu dynamics to the free- 
fermion-like component of the correlated system, in particular since it appears to be quite 
independent on J provided that J < t. Also the free-fermion part seems to exhaust the 
spectra in the overdoped cases > 0.5. Although this observation is not unexpected [9,15], 
the coexistence of spin-fiuctuation and free-fermion timescales at optimum doping has been 
established here for the first time. It is namely harder for other methods, e.g. the high-T 
expansion method, to treat coexisting different timescales. The dual character is a crucial 
property, since the free-fermion part determines to large extent static spin correlations S{q) 
(and charge density correlations N[q))^ interpreted in terms of quasi-Fermi surface [15]. On 
the other hand the low-cu spin dynamics dominates dynamical and static spin susceptibilities, 
x" {q^ io) / id and xl*?)? respectively, hence also the neutron scattering and NMR processes. 

Figure 2 displays dynamical spectra x"((f, c<;)/c<; for fixed = 3/16, but various T and 
nonequivalent q (note that on the 4x4 lattice q = (0, tt) and q = (7r/2, 7r/2) are equivalent). 
In contrast to Fig. 1, high-cu features are suppressed here. Nevertheless, the free-fermion 
part is well separated from the low-cu part only for (f ~ for which one expects a gap in 
the response of free fermions with a well defined Fermi surface. For other q (at given low 
doping) the free-fermion contribution persists at larger > J in the form of a long tail, 
while in the low-cu regime it merges with the spin contribution. The most striking feature of 
Fig. 2 is however strong T dependence of low-u spectra, whereas at the same time the AFM 
correlation length ^ is at most only weakly T dependent. Here ^{T) can be estimated from 
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the static oi" from the ratio of spectra with different q. The above conclusion seems to 
hold for all \q — Q\ < a^~^ ^ where ~ x{Q)- The relevant volume in the (f-space clearly 
increases on doping and exhausts for = 3/f6,4/f6 already the majority of the Brillouin 
zone, while within the same doping regime scaling does not hold for e.g. q = (0,7r/2). The 
variation at < J is x"((f, c<;)/c<; oc f/T, or equivalently S{q,io) is nearly independent on T 
and at < J in the same regime, as also observed in neutron-scattering experiments [f], 
where x"(^) ^ io/T. It should be noted, that the universal scaling x"(^) = fi'^/T) claimed 
by several authors [I] seems to be close to the requirement of S{q^ co) ~ const . 

Spectra discussed above have as a direct consequence T-variation of static ^^^^ 
T < J. We observe pronounced T dependence, e.g. ^ for = 3/16, in a wide 

regime J/3 < T < t for all q within the correlation volume. It should be however noted 
that we are quite restricted in the range of T/J, so that more quantitative conclusions on 
possible power-law (or logarithmic) variation with T are not feasible. 

We can discuss our results in relation to experimental ones obtained in cuprates via the 
NMR and NQR relaxation. The NQR spin-lattice relaxation time Ti and the spin-echo 
decay time T2 for ^"^Cu nuclei are related to electronic spin susceptibilities by [4,3] 



2^ ^ A2t^ x"{q,^o) n 



I 0.69 I 



T2 V 8 g^fil 



1/2 

(4) 



Note that in the evaluation of Ti we neglect the q ^ spin-diffusion contribution, which is 
supposed to be less important [10]. To allow a direct comparison with experiments we choose 
A±[q)^ as proposed in the literature [4] and t = 0.4 eV [7]. Again, J = 0.3 t = 0.12 eV. 

Results for Ti and for the ratio R = T1T/T2 are presented on Figs. 3a, b. For the undoped 
case our results for Ti agree with Ref. [10], but the T-variation for T < J is already influenced 
by finite-size effects (due to sparse density of states at low energies). It is remarkable, that Ti 
appears to be nearly T independent for a broad range of hole concentrations 0.06 < < 0.31 
{ Cfi = 6/16 is not included due to more pronounced finite-size effects caused by a closed- 



shell configuration) . Only for the overdoped systems with > 0.5 the behavior for T < t 
approaches that of a normal Fermi liquid with T^^ oc T. Our results are in agreement, even 
quantitatively without any fitting parameters, with recent remarkable NQR experiments in 
La2-xSrx04 [3], which reveal nearly T and x independent Ti for T > 300A^ and x < 0.15. 
We in fact establish a variation of Ti with doping, which however becomes more pronounced 
only for Ch > 2/16. For optimum doping lower rates T^^ are expected, consistent with the 
data for Y Ba2Cu^07 [2], where for T > Tc Ti is again found to be only weakly dependent 
on T. 

Temperature independent ratio R = T-^T as approximately realized in cuprates, has 
been interpreted as the evidence for the QC behavior of the effective spin system [5]. We 
find quite analogous weak T variation of the ratio within the t — J model, with results 
presented in Fig. 3b. Here undoped c/^ = case is omitted due to inaccurate (finite-size 
dominated) results for x{Q) ^"^^ consequently for T2 obtained on a small system for T < J. 
The origin of R{T) ~ const is however considerably different from the QC scenario, since 
the T2{T) dependence is not connected (in an evident way) with the ^{T) variation. Results 
indicate on a stronger doping dependence, even at low doping. Quantitatively obtained 
values are in reasonable agreement with experimental ones, e.g. R ~ 1700A^ at T = 300A^ 
for YBa2Cu307 {ch ~ 0.23), while R ~ 2400A" for y^aaCMsOe.es [5]. 

In conclusion, we have presented results for the dynamical spin susceptibility obtained 
within the t — J model with the new numerical method, for the first time in the challenging 
regime of low to moderate doping with J < t, as well as T, < J. The most interesting 
finding is the anomalous low-frequency spin dynamics, showing up also in the T dependence 
of static susceptibilities, which is related to NMR relaxation times Ti and as well as to 
neutron scattering experiments. Our results are consistent with nearly T- and cu-independent 
S'((f, cu) (or x"((f, cu) oc io/T) for q around even in the regime where ^ ^ but still 

well in the normal state. Such T dependence is close to the concept of the marginal Fermi 
liquid [6]. 

It is tempting to speculate on the origin of the anomalous low-cu spin dynamics. It is 
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quite plausible to relate it to the dramatic increase of the density of low-energy many-body 
states emerging on doping. The latter shows up in substantially enhanced entropy at low T 
[16], and in the strong carrier scattering as manifested in the resistivity oc T, established 
recently also within the t — J model [14]. An intuitive picture might be that at low T < J 
spin clusters with the characterictic size / ~ ^{T) behave as nearly independent whereby 
their interaction is effectively blocked by doped holes, thus leading to large density of low-cu 
modes. A more coherent theoretical description is clearly missing so far. Our study shows 
that the t — J model still remains a promising starting point for these investigations. 
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Figure captions 

Fig. 1. Dynamical spin-structure lactor S{q = Q^io) at various hole dopings and fixed 
T = 0.2 t < J. Spectra for planar 16-site system are broadened with rj = 0.07 t. 

Fig 2. Dynamical spin susceptibility x"((f, c<;)/c<; at fixed doping = 3/16 for different q 
in the Brillouin zone and various temperatures: T/t = 0.1 (full line), 0.2 (dashed line), 0.3 
(dash-dotted line) and 0.5 (dotted line). 

Fig 3. a) NQR spin-lattice relaxation rate 1/Ti and b) the ratio TiT/T2^ both vs. T/t 
for various dopings c/^. Note that here t = 0.4 eV = 4640 K . 
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